Fractional plateaus of the Coulomb blockade of coupled quantum dots 
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Ground-state properties of a double-large-dot sample connected to a reservoir via a single-mode 
point contact are investigated. When the interdot transmission is perfect and the dots controlled 
by the same dimensionless gate voltage, we find that for any finite backscattering from the barrier 
between the lead and the left dot, the average dot charge exhibits a Coulomb-staircase behavior 
with steps of size e/2 and the capacitance peak period is halved. The interdot electrostatic coupling 
here is weak. For strong tunneling between the left dot and the lead, we report a conspicuous 
intermediate phase in which the fractional plateaus get substantially altered by an increasing slope. 

PACS numbers: 73.23.Hk,72.15.Qm,73.40.Gk 



I. INTRODUCTION 

At low temperature, the charge on an isolated metal- 
lic grain (micronmetric dot) is known to be quantized 
in units of the electron charge e. Even when the grain 
is weakly-coupled to a bulk lead, so that electrons can 
occasionally hop from the lead to the dot and back, the 
grain charge remains to a large extent quantized^. This 
is commonly refered to as the Coulomb blockade?. In the 
opposite limit of perfect transmission between the reser- 
voir and the dot, the average dot charge now depends 
(in a continuous manner) linearly on the applied gate 
voltage and the Coulomb blockade disappears^. How- 
ever, Matveev has shown that a crossover from the linear 
charge-voltage dependence to a Coulomb-staircase func- 
tion occurs for any finite backscattering from the quan- 
tum point contact (QPC) between the grain and the 
lead*. The physics remains qualitatively unchanged by 
increasing the reflection amplitude at the QPC. 

Furthermore, close to the steps, the charge exhibits a 
nonanalytic logarithmic dependence on the voltage due 
to the presence of two spin channels entering the dot, 
resulting in an underlying two-channel Kondo model^. 

Note also that the Coulomb blockade can be smeared 
out by applying an in-plane magnetic fieldA. 

A direct measurement of the average grain charge, 
has been made possible using a single-electron transis- 
tor (SET) which has a sensitivity well below a single 
charge as well as a small input capacitance2i In par- 
ticular, some of the predictions above have been checked 
experimentally and its superiority to conductance mea- 
surements of charge fluctuations demonstrated^. 

Here, we investigate exotic Coulomb staircases with 
fractional plateaus. 

The simplest system we consider comprises two large 
symmetric dots, which can be viewed as an artificial 
molecule, connected to a single reservoir via a single- 
mode QPC (Fig. 1). For a recent review on artificial 
molecules built up with two dots, see Ref. 0. Here, each 
dot is coupled with the same capacitance C g d to a side- 
gate. The term "large dot" implies that the spacing 
A ~ L~ 2 of the energy levels on each dot vanishes com- 
pared to the dot's charging energy E c = e 2 /(2Cs) ~ L~ l 



FIG. 1: A two-dimensional electron gas (2DEG) is coupled to 
two large dots via a single-mode QPC. The number of elec- 
trons Qi on each dot is controlled via the dimensionless gate 
voltage Ni\ The case of interest here is Ni — N2 — N. The 
auxiliary gates can be used to adjust the conductances at the 
QPCs. A SET may probe the average charge on a single dot. 



where Cs w C g d- We already stress that strong tunneling 
between the dots ( "covalent binding" ) is required in order 
to find Coulomb staircases with fractional steps. 

For example, when the interdot transmission is per- 
fect and dot 1 is weakly-coupled to the lead, the inter- 
dot charge fluctuations are so strong that only the total 
charge of the two dots, eQ = e(Qi+Cj2), can be quantized 
(but not the charge eQi on an individual dot). Thus, the 
electrostatic Hamiltonian of the two dots can be rew rit- 
ten more conveniently as (for details, see Refs. Il0lll|) 



H c [N]^^(q~2nY 



2E C N 2 . (1) 



The interdot capacitive coupling is weak in order to max- 
imize the interdot charge fluctuations^. Moreover, the 
symmetric dots are controlled by the same gate voltage 
Vg and N = VaC g d/e. From the electrostatic Hamilto- 
nian, it can be easily inferred that the double dot behaves 
as a single composite conductor of quantized charge 
eQ = 2eQi determined by the total gate voltage 2N. 
When an electron tunnels into the left dot, i.e., Q = 1, 
this implies that a charge e is fluctuating back and forth 
between the dots and clearly Qi=\,2\N\ exhibits steps of 
size 1/2 13 . Moreover, close to a point 2N* = (2n + l)/2 
(n e A/"), the charge states with Q = n (Qi — n/2) 



and Q = n + 1 (Qi = n/2 + 1/2) are degenerate re- 
sulting in (sharp) peaks in the single dot capacitance 
C\ oc dQi/dN (Fig. 2). Similar to the conductance peaks 
for two large dots tunnel- coupled to leadsi2*ii, we then 
observe that strong interdot charge fluctuations produce 
the halving of the capacitance peak period. For an ex- 
perimental proof, see e.g. Ref.hJ 

Based on two-impurity two-channel Kondo models 
(2CKMs) (small dots coupled to leads are described by 
a two- impurity lCKMi£), below we thoroughly analyze 
the evolution of the fractional steps as a function of the 
hopping parameters t\ and t% (Fig. 1). Some aspects 
of the problem will join up with previous works on the 
conductance through a double (large) dot structur e) 1 ?' 11 . 

From here on, we assume that a single orbital channel 
with two spin polarizations a — |,J, enters the double 
dot. Again, we assume that the level spacing on each dot 
(almost) vanishes which means that we consider a contin- 
uous spectrum in each dot and we neglect the mesoscopic 
corrections to the capacitance C g d\ the size of a dot can 
thus exceed the effective Bohr radius (~ (im in Refs. 
I§ IT^) . Temperature will be taken to be zero (T = 0). 



II. WEAK COUPLING WITH LEAD 

Weak tunneling (t\ -C 1) between the lead and the 
composite dot produces corrections to the Coulomb stair- 
case behavior found above. 

More precisely, for perfect interdot transmission {t% — > 
1), we can describe the composite dot in the vicinity of 
the two QPCs by the same field operator ^ ca (x). Ad- 
ditionally, close to a degeneracy point N* = (2n + l)/4, 
only the states with Q = n and Q = n + 1 are allowed 
and, thus, following Ref. 01 the tunneling Hamiltonian 
for this truncated system takes the form: 



Hi 



E 



t 1 *i a {0)V ra (0)S + + h.c. 



(2) 



^ ra stands for the electron operator in the lead, and the 
spin operator S + guarantees that when an electron tun- 
nels into the double dot, the total charge Q only changes 
from n to n + we then have the equalities^ 



= 2Qi = (n + 1/2) 



(3) 



Following the route of the single-dot problem^ 6 , now 
we can identify s~(0) = 5 , J ct (0) v I' rQ (0) as an electron 
pseudo-spin operator acting on the (orbital) indices j = 
r, c and finally recover a 2CKM&±£. The two channels 
are the two spin states of an electron. In particular, 
Eq. Q can be viewed as a local magnetic field hS z with 
h oc (2ra + 1 - 4jV). This results in 



Q 1 „^l±l 0C (2n+l-4iV)ln(|iV- 2 " 



4 ' m 4 

To sum it up, we recover a standard logarithmic form 

2n+ 1 



(4) 
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FIG. 2: Charging energies (+8E) of the "composite" dot as a 
function of N given in units of E c ; ti is small. Each eigenstate 
with Q = n gives rise to a parabola. The solid lines corre- 
spond to V2 — and dashed lines to increasing couplings. 
For T\ = 1 — ti — > 1, Q is quantized and for symmetric dots 
this guarantees Qi = n/2 until ri — > 1 (Eqs. 111U24I 1. 

for the capacitance peaks. 

We now discuss the situation in which the interdot 
tunneling is strongly decreased, (iijia) -C 1. Each dot 
is described by its own operator and the Coulomb 
term should be written in a more common way asiii^: 

H\[N\ + H 2 C [N] =E C Y] (Qi - N) 2 - 2E C N 2 . (6) 



When t% — > 0, we converge to a single-dot problemii&: 
Qi is quantized and we could not use Qi ~ Q/2 in Eq.Q! 
The degeneracy points now occur for N* = (2k + l)/2 
(k € Af) and obviously the period of the capacitance 
peaks then doubles (Fig. 2). 

As soon as ti is finite (t^ ~ t\) and N w 1/2, we 
propose to modify the tunneling Hamiltonian as: 

H t =J2 (*l*la(0)*ra(0) + h^jL^ la {L) + h.C.) (7) 



hs-(0)S+ + t 2 s-{L)S+ + h.c 



where S±(Sz) emphasizes that the charge on dot 1(2 ) 
only changes from to 1; For more details, see notellSl 
Again, the index j = 1,2, r -which designates the loca- 
tion of an electron in the setup- in the "i>j a operator can 
mimic an internal "orbital" degree of freedom. It is then 
straightforward to define two spin operators at x = (0, L) 
acting on the orbital space, similar as in Ref.Q: 



«a(0) = *L(0)*™(0) 



(8) 



SCt = Ci - Cga <x -\b.{\N 



(5) 



This two-impurity (two-channel) Kondo model is partic- 
ularly convenient to revisit the behavior of charge fluc- 
tuations close to the degeneracy points N*; the crucial 



3 



point being that a finite bare coupling t 2 (like ii) will be 
strongly renormalized at low temperatures^. 

At the fixed point (T — 0) and, e.g., close to the degen- 
eracy point iV* = 1/2, the two dots will merge into one 
and therefore by analogy to Eq. @ we must correctly 
reidentify 



= (j + l/2) + S z = 2Q 1 , 



(9) 



where j = (0; Moreover, the Coulomb term in the 
fixed-point basis takes the form hS z where h cx (1 — 2iV=F 
2kT 2 ) for j = (0; 1); 7a = (t 2 ) 2 and k > 0. Away from the 
point N* — 1/2, second order perturbation theory in t 2 
is accurate, and we have taken into account the relative 
energy shift between even and odd Q-statesifliii: 



SE oc -4% In 2. 



(10) 



Similar to Eq. 10} , we are thus led to (for j — 0, 1 respec- 
tively) 



b(N 



kT 2 
kT 2 



N ■ 
N - 



kT 2 
kT 2 



(11) 



b > is a parameter which is inversely proportional to the 
Kondo energy scale. By continuity, a tiny step appears 
at Qi = 1/2, and the single-dot capacitance peaks are 
already split by ~ 2kT 2 (Fig. 2). 

The progressive pairing of the capacitance sub-peaks 
close to t 2 = 1 will be studied later (Eq. (24)). 



III. STRONG COUPLING WITH LEAD 

Now, we mainly consider the case where all the junc- 
tions have conductances close to 2e 2 /h, i.e., reflection 
amplitudes are small {r\\r 2 ) <C 1. 

In this case, the whole system can be viewed as a 
single conductor and, for convenience, we will use the 
unique field operator 1 $ ra (x)^. We can write ^ ra (x) — 
CTq>(ikpx)ty +a (x) + exp(— ikpx)^ - a {x), <5 ' + a and 
describe right- and left-moving fermions respectively. 
The kinetic energy obeys 



Hi. 



ivf 



-2L 



dx (vl a d x * +a - vl a d x y- a ), (12) 



Vf being the Fermi velocity. The backscattering term(s) 
takes the standard form: 



(13) 



and interactions in a grain are embodied via the general 
Coulomb Hamiltonians H^Nx] + H C [N 2 ], in Eq. ©. 

At low energy, we proceed with this model by bosoniza- 
tion of the one-dimensional Fermi fields 6 . In those vari- 
ables, the kinetic energy yields a separation of the spin 



and charge and the resulting Hamiltonians have plasmon- 
like excitations. Here, d x <j)j with j = (c, s) measures fluc- 
tuations of charge/spin density and Uj = d x 6j being its 
conjugate momentum. The Coulomb Hamiltonians take 
the forms (We could also employ Eq. Q): 

H l cWi] = ^ (&(0) -&(£)- ^/f- 



H. 



w - 2 M 



C (L)-0 C (2L)- A /-7V 2 



V, ) -E e Ni*(U) 

2 

- E r N 2 2 . 



To minimize the Coulomb energies when the transmis- 
sions at the two QPCs are both perfect, we easily recover 
that the dot's charges evolve continuously (linearly) as a 
function of the gate voltages^: 



! = —(&(0) =N! 



(15) 



Q-2 



V2 



(<j> c {L) - <f> c (2L)) = N 2 



Remember that for r\ — r 2 = 0, the Coulomb blockade 
physics is totally suppressed. In our geometry there are 
no charge fluctuations at x = 2L and then <fi c (2L) = est. 

Furthermore, following the traditional route of the 
single-dot problem for this regime 4 * 6 -, the backscattering 
term may be rewritten as 



H h = 



y/^aE c v F 



y/jaE c v F 



it a 



An cos^A^ + N 2 )) cos (VW> s (0)) Tlx (16) 



4r 2 cos (ttN 2 ) cos V2tt0 s (L) T 2x 



Since the charge fluctuations on each dot cannot depend 
on the precise size of a dot, we must equate <f> c (2L) = 
2k F L/V2~w and rescale cj> c (0) -> <j) c (0) + 2k F L/V2~K. 
Here 7 obeys 7 — e c where C rs 0.5772... is the Eulcr- 
Mascheroni constant and a is a short-distance cutoff. We 
have introduced two commuting impurity spins T\ and T 2 
(which here are not related to the charge on each dot). 
Clearly, the T\ x and T 2x spin-operators both commute 
with the Hamiltonian and must be simply identified as c- 
numbers, i.e., T\ x = 1/2 (or —1/2) and similarly for T 2x . 
Eq. I|16l) must be viewed as an extension of the 2CKM 
at the Emery-Kivelson lign22. 

To compute the correction to the average dot 
charge(s), here we must "de-bosonize" the problem as 4,6 



+ 



/JlK =(V»(o) + ^(o))Ci 
= (V;(L) + V t (i))C2. 



(17) 



£1 and £2 are two Majorana fermions, and the Kondo ex- 
changes above read J\ x 
and J 2x 



Ar x ^a^E c v F cos(7r(A^i + N 2 )) 

4r 2 yJa^E c v F cos(ttN 2 ). In the absence of an 



applied magnetic field, there is no net magnetization and 
no spin current on the whole region [-L; L] and, thus, we 
have approximated^ ip(L) ~ exp(i\/27r^i s (i)) (For more 



4 




N ]= N 2 



FIG. 3: Dot's differential capacitance for small n = 0.4. For 
T2 ~ the system behaves as a single composite, conductor 
whereas for r2 — ► 1 we must recover a single- dot problem. 
For f2 ~ ri < 1 the system cannot decide between those 
two ground states giving a "3-peak" capacitance profile, i.e., 
unstable fractional steps (Fig. 4 and Eq. 1221 1: Charge fluc- 
tuations are important at N = 1/4 and at N = 1/2 as well. 
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FIG. 4: Evolution of the fractional plateaus for small n. 
(Note e.g. that for r2 = we have only taken into account 
the main leading term when taking the derivate of Eq. (19), 
which explains the "slightly" negative slope in the middle of a 
plateau). The solid line is for n = 0.2 and r2 = 0, the dashed 
line for r\ = r-z = 0.15 (fractional plateaus now acquire a 
positive slope) and the dotted line for n = 0.3 and r2 — » 1. 



B. 



explanation, see Ref. [fj). The fermionic model here gen- 
erates two Kondo resonances 



Jl X _ ^l {2ri) 2 cos 2 (7T{Ni+N2)) (lg) 



Airavp 

2 



T 2 = -^ = ^I(2r 2 ) 2 cos 2 UiV 2 ), 

and all the quantities of interest will be now inferred from 
the quantum correction of the ground state energy 

6E = -^]n(E c /T 1 )-^lii(E c /T 2 ), (19) 

which implies that impurities are independently screened. 
Let us discuss the case of symmetric dots: Aq = N 2 = N. 
The correction to the average charge on each dot 8Qi 
and the dot's differential capacitance 8Ci obey: 5Qi — 
Qi- N cx —d5E/(E c dN) and 5C t cx dSQi/dN. For 
clarity's sake, results have been summarized in Figs. 3,4- 



r-2 



For a double dot connected by a reflectionless constric- 
tion <— r 2 <C ?'i( < C 1), using the formulas above, we 
easily recover fractional charge plateaus with steps 1/2 
and capacitance peaks with halved period. Again, the 
double dot behaves as a single composite conductor of 
quantized charge Q w 2Q±. In particular, we predict that 
the logarithmic sing ularity 5C l cx -h\(\N - \\) should 
be observed at any value of t\ ^ 1 (as nicely illustrated 
in Figs. 3 and 4). 



In the opposite limit rq <tC r% — ► 1, the interdot con- 
striction considerably impedes the charge spreading be- 
tween the dots. Qi becomes an integer-valued operator 
describing electrons which tunnel into the dot 2 (Eq.©) 
and charge fluctuations in dot 1 closely resemble the ones 
of a single dot which is strongly-coupled to one lead: 

H b cx ri(-l) Q2 cos(ttA0cos (V2n(f> s (oj\ T lx . (20) 

The Kondo energy scale 



Fl = ^2(2r 1 ) 2 cos 2 ( 7 rA^), 

7T 



(21) 
and 



is identical to the one of the single-dot problem 4 
assuming rq 7^ 0, Q\ becomes also quantized. 

The small term t 2 ^l a (L)^ la (L)S2 + h.c. here mostly 
produces slight charge fluctuations in the dot 2, and 
6C 2 cx -ln(|iV- 1/2Q. 



C. n « r 2 < 1 

For rq w r 2 , a strong opposition between the single- 
dot (Qi is quantized for r 2 ^> r%) and the composite- 
dot ground-state (Qi = Q/2 for r\ ^S> r 2 ) arises giving a 
fascinating "hybrid" regime where the fractional plateaus 
become gradually destroyed by acquiring a positive slope 
(Fig. 4); Close to N = 1/2, exploiting Eqs. ijTSj) and l {Tg| l. 
we can approximate (f(TZ±) = ln7?.i + const.) 

SQt cx (JV - 1/2) (^1/(^1) - ^ 2 ln(|AT - 1/2|)), (22) 

then inducing an exotic "3-peak" capacitance profile; 
Hi = ( r i) 2 (inset in Fig. 4). The central peak becomes 
more pronounced by slightly increasing r 2 , whereas the 
external peaks only depends on rq (as long as r 2 <C 1). 
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D. n — ► 1 and r2 1 

It is worthwhile to compare with the case r\ — » 1 
and r 2 <C 1. Here, Q = y^2/n(j) c (0) — n must be 
an integer-valued operator which guarantees Qi = n/2. 
The fractional plateaus remain by decreasing the inter- 
dot coupling and only their widths progressively reduce: 
N*(n = 1)-N*(n = 0) = l/2-2r)H 2 ln(l/K 2 j\ 77 > is 
a constant parameter. More precisely, for Ni = A 2 = N , 
it is easy to rewrite the backscattering term as: 



F 6 = ^^4r 2 cos(^)cos(^(L))T 2a , (23) 

to Eqs. CQ,© which then produces a Kondo energy scale 
r 2 = E c j(2r 2 ) 2 cos 2 ( ) /tt, and then a relative energy 
shift 5E oc IZ2 ln(l/7?. 2 ) between even and odd statesiSaii. 
This engenders that the positions of the capacitance (sub- 
)peaks (furnished by Eq. 1) are shifted as 

N* = (2n + l)/4 + (-1)"^2 ln(l/7e 2 ). (24) 

The capacitance (sub-)peaks are not equally spaced any- 
more and progressively pair around the points N* — 
(2n + l)/2 (Fig. 2). Finally, we have checked that integer 
plateaus become more prominent: N*(n — 2) — N*{n — 
1) = l/2 + 2^ 2 m(l/ft 2 ). 

IV. CONCLUSION 

In closing, based on two-impurity two-channel Kondo 
models, we have presented a detailed discussion on the 



evolution of the fractional plateaus as a function of the 
hopping parameters t\ and t 2 for a double-dot coupled 
via a single-mode QPC to a reservoir. 

Again, for perfect interdot transmission, Coulomb steps 
of size 1/2 occur for any finite backscattering between the 
lead and the left dot. When an electron enters the arti- 
ficial molecule, a charge 1 is fluctuating back and forth 
between the two dots. We are hopeful that this can be 
observed via capacitance measurement s" 2 4ff 

Substantially decreasing the interdot coupling in- 
evitably restores the single-dot Coulomb blockade and 
the capacitance peak period doubles. 

For strong coupling between the lead and the left dot 
(ri -C 1), we find a striking intermediate range (r 2 ~ ^1) 
where the fractional steps become progressively unsta- 
ble, i.e., show an increasing positive slope; this happens 
due to the strong competition between a single-dot and 
a composite- dot ground state. On the contrary, when 
r% — > 1, Q must be quantized and Qi = Q/2; the frac- 
tional steps persist. 

For asymmetric dots, e.g., with different gate-dot ca- 
pacitances, we report that the Coulomb staircase with 
halved steps is gradually altered. 
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